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Abstract. Anomalous dynamics characterized by non-Gaussian probability distribu¬ 
tions (PDFs) and/or temporal long-range correlations can cause subtle modifications of 
conventional fluctuation relations. As prototypes we study three variants of a generic 
time-fractional Fokker-Planck equation with constant force. Type A generates su¬ 
perdiffusion, type B subdiffusion and type C both super- and subdiffusion depending 
on parameter variation. Furthermore type C obeys a fluctuation-dissipation relation 
whereas A and B do not. We calculate analytically the position PDFs for all three 
cases and explore numerically their strongly non-Gaussian shapes. While for type C we 
obtain the conventional transient work fluctuation relation, type A and type B both 
yield deviations by featuring a coefficient that depends on time and by a nonlinear 
dependence on the work. We discuss possible applications of these types of dynamics 
and fluctuation relations to experiments. 
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1. Introduction 

Understanding fluctuations far from equilibrium defines a key topic of nonequilibrium 
statistical physics. A new line of activities started about three decades ago by discovering 
different forms of fluctuation relations (FRs) which generalize fundamental laws of 
thermodynamics to small systems in nonequilibrium; see Refs. [1, 2, 3, 4, 5, 6, 7, 8] 
for reviews and further references therein. More recently these laws got unified by 
over-arching schemes, most notably the deterministic dissipation function approach by 
Evans and coworkers [1], and by stochastic thermodynamics [9, 10, 7, 11]. The latter 
theory starts from defining entropy production on the level of individual trajectories 
in stochastic models such as Langevin and master equations. Given that stochastic 
thermodynamics is based on rather simple Markov models one may ask to which extent 
FRs derived from it are reproduced if the dynamics is more complicated. In our paper 
we address this problem by testing FRs for stochastic dynamics that is anomalous due 
to non-Markovian dynamical correlations and/or strongly non-Gaussian PDFs. 

Anomalous dynamics has been observed in many experiments and is widely studied 
by the theory of anomalous stochastic processes [12, 13, 14, 15, 16, 17]. A characteristic 
property of anomalous dynamics is that the mean square displacement (MSD) grows 
nonlinearly in time yielding anomalous diffusion in the long time limit [15]. In contrast, 
Markovian dynamics like Brownian motion generates a MSD that increases linearly 
for long times. If the MSD grows faster than linear one speaks of superdiffusion, 
if it grows slower than linear one obtains subdiffusion. There are many different 
ways to model anomalous stochastic dynamics such as continuous time random walks 
(CTRW) [18, 19, 20, 12], generalized Langevin equations [21, 13, 22, 23], Levy flights 
and walks [24, 17], fractional diffusion equations [16], scaled Brownian motion [25, 26] 
and heterogeneous diffusion processes [27], to name a few. 

The study of FRs for anomalous stochastic processes appears to be rather at 
the beginning: Crooks and Jarzynski work relations as well as transient and steady 
state fluctuation theorems have been confirmed for non-Markovian Gaussian dynamics 
modelled by generalized Langevin equations with memory kernels, given specific 
conditions are fulfilled [28, 29, 30, 31]. These results have been reproduced and 
generalized by a stochastic thermodynamics approach [32]. For non-Gaussian PDFs 
generated by Langevin equations with non-Gaussian noise, such as Levy noise or 
Poissonian shot noise, violations of conventional steady state and transient FRs have 
been reported [33, 34, 35, 36, 37, 38]. For a CTRW model with a power law waiting 
time distribution it was found that the steady state FRs may or may not hold depending 
on the exponent of the waiting time distribution [39]. Computer simulations of glassy 
dynamics exhibiting anomalous diffusion also showed violations of transient fluctuation 
relations [40, 41], In [42, 43, 44] several of the above types of stochastic dynamics 
including fractional Fokker-Planck equations were considered. It was found that the 
validity of fluctuation-dissipation relations [45] for a given anomalous stochastic process 
plays a crucial role for the validity or violation of conventional FRs. 
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In this article we test transient fluctuation relations (TFRs) for a class of anomalous 
stochastic processes that so far has not been in the focus of investigations, which are 
time-fractional Fokker-Planck equations (FFPEs). Such equations model the emergence 
of non-Gaussian PDFs by using power law memory kernels via time-fractional derivatives 
[46]. They need to be distinguished from equations modeling correlations in space 
via space-fractional derivatives as they naturally arise, e.g., for generating Levy flights 
[12, 42], FFPE can be derived from stochastic equations of motion either by CTRWs 
[12, 16] or by subordinated Langevin dynamics [47]. Quite a variety of them have been 
studied in the literature, both from a purely theoretical point of view and with respect 
to applications to experiments: Prominent examples are fractional Klein-Kramers 
equations that were used to analyse biological cell migration data [48, 49, 50]. Another 
type was designed to model the dynamics of tracer particles in random environments 
[51]. Closely related time-fractional diffusion equations [21, 52, 12] have been used to 
model a variety of different processes, from diffusion in crowded cellular environments 
[53, 15] to geophysical and environmental systems [14]. They have also been derived for 
weakly chaotic dynamical systems [54, 55]. A bifractional diffusion equation famously 
reproduced the spreading of dollar bills in the United States [56]. 

Our paper is structured as follows: In Section 2 we discuss three types of FFPEs 
which differ from each other in terms of their anomalous diffusive properties, and by 
whether or not they fulfill fluctuation-dissipation relations. We solve these models 
for their position PDFs and study their properties both analytically and numerically. 
In Section 3 we test the (work) TFR for our three models by analytical asymptotic 
expansions and by numerically plotting the results. We conclude with a summary and 
an outlook towards physical applications in Section 4. 

2. Time-fractional Fokker-Planck equations 

This section introduces to three different types of FFPEs: We first outline how starting 
from stochastic dynamics a FFPE generating superdiffusion can be constructed in 
the form of an overdamped Langevin equation with correlated noise. Our argument 
illustrates how a time-fractional derivative naturally emerges from modelling power law 
time correlation decay. The other two types of FFPEs that we consider have already 
been derived in the literature from CTRW theory and are either subdiffusive or exhibit 
a transition from sub- to superdiffusion under parameter variation. We analytically 
calculate the first and second moments for all three models, which enables us to check 
for the validity of the fluctuation-dissipation relation of the first kind (FDR1). We also 
comment on the Galilean invariance of our models. We then analytically calculate the 
position PDFs of all FFPEs and study the solutions numerically by plotting the results. 
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The study of an overdamped Langevin equation for the position x(t) of a particle on 
the line driven by a correlated stochastic process and an external force allows to gain 
insight into the origin of a superdiffusive FFPE. Our Langevin equation of interest is 
given by 


dx 

dt 


Fp 

m7 Q 


+ v(t), 


( 1 ) 


where F 0 denotes a constant external force, a friction coefficient and m the mass of 
the particle. We assume that v(t) is a stationary correlated stochastic process with zero 
mean (v(t)) = 0 and a power-law correlation function 


(v(t) u(ti)) 


Kg 1 

r(a-l) 


( 2 ) 


with f < a < 2, gamma function F and generalized diffusion coefficient K a . Note that 
we do not further specify the noise. Following the pseudo-Liouvillc hybrid approach of 
Balescu [57, 58] (see Appendix A) one obtains the following exact result in Eq. A.8 for 
the PDF ffx , t)\ 


(jf t + A) f(x, t) = ^~2 J dt 1 (v(t)v(t i) / (x - A(t, h), h) ) (3) 

0 


t 

with vq = F 0 /(m'y a ) and A(t,ti) = vq( t — ti) + / dt 2 v(tf). This exact equation is 

ti 

non-local in time (i.e. non-Markovian) and non-local in space. We now make a local- 
in-space approximation by neglecting the term of the fluctuating displacement A (t,ti) 
on the right hand side of the probability density /. Such an approximation seems to 
be reasonable in the long time and large space asymptotic limit if the drift and velocity 
fluctuations are weak enough. This assumption results in the following non-Markovian 
Fokker-Planck equation: 

t 

+ f( ' Xl ^ = ^ ^ ^ ^ 

0 

Insertion of the correlation function of velocities Eq. 2 into Eq. 4 leads to 


(d d\ 

\dt +v,> dx) 


f(z,t) 


d 2 Kg 

dx 2 r(a — 1) 


t 

j dp (t - ti ) a ~ 2 f(x,ti) . 
o 


( 5 ) 


The integral on the right hand side matches to the Riemann-Liouville (RL) fractional 
integral of order p given by [59] 


Jt 9(t) = D t M g{t) 


1 

foo 


dr (t - tY 1 g{r) 


o 


( 6 ) 
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with p > 0 and fi = a— 1 for Eq. 5. We also introduce the definition of the RL fractional 
derivative of positive order 

d n 

D t 9(t) = ^ J?- 11 g(t ) (7) 

with fi > 0 , n — [/i] + 1, where [...] refers to the integer part of the given number. 
Applying Eq. 6 to Eq. 5 gives us our Erst type of FFPE that we denote as 


type A: 


df A (x,t ) 
dt 


d_ 

dx 


v 0 - K a D 


l—a 
t 


dj 

dx 


f A (x,t) , 1 < a < 2 . 


( 8 ) 


To show the relation of this equation with previous works we put Vq = 0. Then it can 
be written as 

d 2 d 2 

Qpf( x i *) = K « D t a f( x i t) ,1 < a <2. (9) 

This equation was called a fractional wave equation in the seminal paper of Schneider 
and Wyss [52] and has also been derived for a long-range correlated dichotomous 
stochastic process [60] from a fractional Klcin-Kramers equation [48] and from a 
generalised Chapman-Kolmogorov equation [61]. The solution of this equation has 
been studied in detail in [62] where it was called a fractional kinetic equation for sub- 
ballistic superdiffusion. The equivalent form of this equation using the Caputo fractional 
derivative was investigated in [63]. 

Our presentation above illustrates how a FFPE can be derived from a Langevin 
equation with power-law decay in the velocity correlation function. It furthermore 
demonstrates that a fractional derivative provides the natural mathematical formulation 
to model equations containing power law memory kernels. 


2.2. Definition and properties of fractional Fokker-Planck equations 


In addition to type A FFPE Eq. 8 we consider two further types of FFPEs. Both 
have been derived from CTRW theory [18, 19, 20, 12]. Note that the underlying 
stochastic dynamics and the derivation of these two FFPEs are very different from what 
we presented for type A above. Indeed, both type B and type C are essentially (almost) 
Markovian models, in contrast to type A. Our two new FFPEs describe subdiffusion 
under the influence of a constant external force and naturally appear in physical systems 
where diffusion is slowed down by deep traps [12, 20, 64], The difference between these 
two types arises from the position of the fractional RL derivative with respect to the 
diffusive and drift part of the equations and the range of the anomaly parameter a. Our 
second FFPE is defined as 


type B: 


df B (x,t ) 
dt 


d_ 

dx 


Vo - K a D]~ a 


dj 

dx 




( 10 ) 


For type C FFPE the RL fractional derivative is also included in the drift term: 


type C: 


dfc(x,t ) 

dt 


d_ 

dx 


A a v 0 D}~ a 


K D 1_a — 
KaVt dx 


fc(x,t), 


( 11 ) 
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where A a has a dimension of time to the power of 1 — a. Note that type B and type C 
FFPEs are defined for 0 < a < 1 whereas for type A FFPE a is in the range 1 < a < 2. 
For all three FFPEs we use the initial condition fA,B,c(x,t = 0) = h(a:). By means of 
Fourier and Laplace transforms 


f(k) = I dxe ikx f(x\ f(s) = I dte~ st f(t) 

—oo 0 

a solution of Eqs. 8, 10 and 11 can be obtained in Fourier-Laplace space as 

/n,s( ,s) g Vo _|_ K a k 2 s l ~ a ’ 


( 12 ) 


(13) 


fc(k,s) = 


(14) 


s + Araiks 1 - 01 + K a k 2 s l ~ a ’ 
where the fractional derivative Dl~ a f(t ) transforms to s 1 ~ a f(s). The solutions of type 
A and type B FFPE only differ in the range of a as defined above. The representation in 
Fourier-Laplace space allows the calculation of moments by differentiation with respect 
to k : 

d n f(k,s ) 


(x n (t)) = /T 1 (*) 


dk r 


(15) 


k—0 ■ 


After Laplace inversion one obtains the first two moments and the central second 
moment for Sx = x — (x) of type C FFPE defined in Eq. 11 [20] 

A a v 0 t a 


(x)c 

(A 

<(^) 2 >C = 


r(a + i)’ 

2IiA c 


>c = 


+ 


2 A 2 y 0 t 2a 


r(a + 1) T(2a + 1) ’ 


2K a t c 


+ Alvlf 2a 


(16) 

(17) 

(18) 


r(a + l) “ u L r ( 2 « +!) T(a + 1 ) 2 . 

These results show that the FDR1 [45, 43] (x(t)) c ~ (x 2 (f))^ 0=0 is valid for 
type C. Interestingly the external force influences the second central moment ~ 
i>Qt 2a . Technically this is due to the coupling term Uo iks l ~ a in the Laplace-Fourier 
representation of Eq. 14. The first moment increases sublinearly despite the constant 
external force. This can be interpreted as a partial sticking effect of particles [65]. By 
contrast, the second central moment shows a crossover from ~ t a to ~ t 2a . Thus, for 
r>o 7 ^ 0 type C switches from a subdiffusive behavior of the second central moment for 
0<a<l/2toa superdiffusive behavior for 1/2 < a < 1 [12]. 

Analogously, the moments of type A and type B FFPEs of Eq. 8 and Eq. 10 are 
obtained as [20] 


(x)a,b 
(■ x 2 )a,b 

((Sx ) 2 )a,b 


Vo t, 


2Kj a 
r(a + 1 ) 
2 K a t a 
r(cr + 1) 


+ v 2 0 t\ 


(19) 

( 20 ) 


( 21 ) 
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In both cases the first moment only depends on vq and increases linearly in time. The 
second central moment shows a superdiffusive and subdiffusive increase ~ t a for type 
A and type B FFPE, respectively. In contrast to type C FFPE, the second moment of 
type A and type B FFPEs is without any coupling to no- In addition, type A and type 
B FFPEs break FDR1 between the first ( x(t)) a,b and the second moment (x 2 (t)) v ^°. 
In both cases this is what one should expect according to the definition of both models: 
Type A is based on the Langevin equation 1 where the fluctuation-dissipation relation 
of the second kind (FDR2) is broken by construction. Note that FDR2 establishes 
a relation between the noise and the friction [45]. The breaking of FDR2 suggests a 
breaking of FDR1 as was shown for Gaussian stochastic processes in [43]. For type B 
the fractional derivative acts only on the diffusion term in Eq. 10 thus breaking FDR1 
while for type C it acts simultaneously on both the drift and the diffusion terms in 
Eq. 11 hence preserving FDR1. 

A second difference between these FFPEs consists in their behavior under Galilean 
transformation. With X = x — Vq t and T = t the PDF /(x, t ) is transformed to D(A", T). 
The coupling of the fractional RL derivative to the vq drift term of type C FFPE in 
Eq. 11 breaks Galilean invariance. However, type A and B FFPE of Eq. 8 and Eq. 10 
fulfill Galilean invariance in the long time and large space limit [20, 12, 66], where they 
can be written as 

8SIa,b(X,T) t , ^_ a 8 2 n A , B (X,T) 

dT ~ ° T dX 2 ’ 1 

This means that in this limit breaking or preserving FDR1 corresponds to preserving 
respectively breaking Galilean invariance in the case of these FFPEs. This property will 
be exploited in the next subsection where we discuss analytical and numerical solutions 
of our three types of FFPEs. 


2.3. Analytical solution of time-fractional Fokker-Planck equations 


Type C FFPE: Fourier inversion [42] leads to the solution of type C FFPE in (x , s) 
space: 


0 a —1 


fc(x,s) = 


^A 2 Q v 2 0 + 4K Q s< 


exp 


A a v 0 x 
2 AT 


— \x 


■sfAlvl + AK a 
2 AT 


(23) 


In this case, a solution in (x, t ) space can be given as a superposition of the a = 1 
Gaussian solution with a Levy kernel [67, 12], However, for numerical analysis we apply 
a direct numerical Laplace inversion of Eq. 23. 


Type A and B FFPE: Analogously to Eq. 23 the solutions of type A and type B FFPEs 
can be calculated in (x, s) space with A a v o —> vq s q_1 to 


~,a— 1 


Ia,b(x,s) = 


\J XqS 2 " -2 + 4AT q 


: exp 


cx .—1 

Vq S X 


2 AT 


— x 


7 


v^s 2a ~ 2 + 4 AT 


2 AT 


(24) 
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As the FFPEs of type A and type B are Galilean invariant in the long time and large 
space limit, the solution for vq = 0 allows the exact calculation of the PDFs with drift v 0 
in this limit [20, 12], which becomes approximate otherwise [66]. The solution to Eq. 22 
is well-known [12] and is given using a Fox //-function (see Appendix B for definitions). 
Thus, applying Galilean transformation and replacing x with x — u 0 t gives solutions of 
type A and type B FFPEs in ( x , t ) space as 


1a,b(x, t) 


1 H U} 

l 

IT 

i 

o 

C-F 

(1- a/2, a/2) 

11 

V K a t a 

(0.1) 


(25) 


These approximate solutions in terms of shifted Fox functions are the basis for our 
further analysis of type A and B FFPEs. 


2-4- Numerical analysis of time-fractional Fokker-Planck equations 

Numerical methods are required to study the analytical results given in form of Fox 
//-functions of type A and type B FFPE and in Laplace space for type C FFPE. 


Type A and type B FFPE: The series expansion of the solution /a, s ((£,/) of Eq. 25 as 
given by Eq. B.3 is used for numerical evaluations, 


fA,B{x,t) 


1 

V4 


OO 


£ 


(~iy 

j\T{l-a(j + l)/2) 


f {x-v 0 t)' 2 

V K a V* 


i/2 


(26) 


with 1 < a < 2 for type A FFPE and 0 < a < 1 for type B FFPE. The series is 
evaluated with multiple-precision arithmetic. 


Type C FFPE: Direct numerical Laplace inversion is applied to Eq. 23 to obtain the 
probability density function fc{x, t). Here we use a multiple-precision algorithm for the 
Laplace inversion based on Talbot’s method [68, 69]. 

Typical behavior in space and time: Fig. 1 shows the time development of the solutions 
f(x, t ) of the three FFPE types for different times t — 1,2,4,8. Parameters were selected 
as A a v o = 1 and K a = 1, the anomaly index a was chosen from a G (0.4, 0.6...1.6). The 
first row shows the Gaussian limit a —>■ 1 for all three types. In this normal diffusive case 
the PDF is spreading with \j2K\t and its center is moving according to vpt. The PDFs of 
type A (left column) and type B FFPE (middle column) preserve this constant drift for 
a^l. However, the shapes of the PDFs of both models immediately change profoundly 
showing characteristically different types of non-Gaussian behavior: For type A the 
PDFs spread superdiffusively with the variance of Eq. 21 by exhibiting a double-peaked 
structure with a dip in the middle. Qualitatively, the highly characteristic double¬ 
peak structure is explained in [58]: The propagator of type A decays asymptotically 
faster than the Gaussian, cf. Eq. B.5. However, since two maxima move away from the 
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origin in the opposite directions, superdiffusion is possible in spite of the thin tail of 
the propagator; see also Eq. 9 [62], Note that there are cusp singularities in all three 
models for a / 1, in contrast to the smooth behavior of the Gaussian PDF shown in 
the top row. In the Galilean invariant cases A and B the propagators are symmetric 
with respect to their cusps, which are translated with velocity no = 1, as it should be. 
For the Galilean non-invariant model C the propagator is asymmetric with respect to 
its cusp, which stays fixed at the origin [12], 


3. Work fluctuation relations for fractional Fokker-Planck equations 


3.1. Definition of fluctuation relations 

Using the results of the previous section, we now study the probability distribution 
p(W, t) of the mechanical work W = —F 0 x generated by the constant external field 
F 0 . For a constant field the probability distribution of work p(W,t) is related to the 
probability distribution f(x,t) of positions by the simple scaling transformation 


p(W,t) 



(27) 


It is the main aim of this work to study the TFR of the work PDFs defined by the 
logarithmic fluctuation ratio 


R (W,t) := log 


p(W,t) 
p(~W,t ) 


(28) 


for the three types of FFPEs. All three FFPE types reduce to a normal Gaussian process 
with drift for a —)• 1. For a Gaussian PDF the ratio R is trivially given by the ratio of 
the first and second central moment, i.e. 2(W)/(5W 2 ) [43]. Thus one obtains a normal 
or conventional fluctuation relation for a = 1, 

R(W,t) |«=i = W = W/{k B T). (29) 

Co-fll 


with a linear increase in W that is independent of time as it has been found for a large 
class of systems [1, 2, 3, 6, 7, 8]. The last expression has been obtained by using the 
Einstein relation K a = k B T / (m^a) with temperature T, Boltzmann constant k B and the 
definition v 0 = F 0 /( ?7 Ga)- The general case for a 1 is studied in the next subsection. 


3.2. Fluctuation relations for fractional Fokker-Planck equations 

Type C FFPE: For this type the fluctuation ratio can be studied analytically [42]. 
With Eq. 23 R(W,f) is given in Laplace space by 


Pc{W, s ) 
Pc(~W,s) 


fc(W/F 0 ,s) 
fc(-W/F 0 ,s ) 


exp 


A a v o 



F 0 K a 


(30) 







p(x,t) p(x,t) p(x,t) p(x,t) 
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Figure 1. Time development of PDFs for type A FFPE (left column), type B FFPE 
(middle column) and type C FFPE (right column) for different values of a (rows) and 
time points t = 1,2,4, 8. Parameters were selected as K a = 1, vq = 1 and A a v o = 1. 
Whereas superdiffusive type A FFPE (left column) and subdiffusive type B FFPE 
(middle column) show a drift and spreading of the PDFs with typical non-Gaussian 
structures for a / 1, type C FFPE (right column) displays a spreading of the PDFs 
together with stickiness to the origin. 


As the right side is independent of the Laplace variable s, the Laplace inverse of the 
PDFs can be calculated directly after multiplication with pc(—W, s ). Thus, despite the 
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Pc(W,t) 
Pc(~W,t ) 


AgV 0 
F 0 K a 


W. 


(31) 


This result based on the Laplace transformed ratio of pc(W,s) seems to be surprising 
with respect to the complex form of the PDF in Laplace space and the asymmetric 
sticking behavior at the origin of the PDFs as illustrated in the right column of Fig. 1. 
The right side of Fig. 2 shows the numerical calculation of the fluctuation ratio which 
is linear and constant for all times in agreement with the given analytical result. 

We remark that a normal TFR for type C can also be obtained with the use of the 
subordination principle: Indeed, it is known that the fractional kinetic equation C can 
be derived from the coupled Langevin equations for the motion of a particle [47, 70, 42] 


dx(u) 

du 


— + ««) 
my 


dt(u) 

du 


t(u) 


(32) 


where the random walk x(t) is parameterized by the random variable u. The random 
process £(u) is a white Gaussian noise, (£(w)) = 0 , {£(u)£(u')) = 2 k b TS(u — u')/(rn , y), 
and r{u ) is a white stable Levy noise, which takes positive values only and obeys a 
totally skewed a-stable Levy distribution with 0 < a < 1. The PDF f(x,t ) of the 
process x(t) is then given by 




dufi(x, u)h(u, t ) 


(33) 


where fi(x,u) is a shifted Gaussian PDF with drift, and h(u,t) is the inverse one-sided 
Levy stable density [67]. It is then easy to show that the linear normal TFR Eq. 31 
holds due to Gaussianity of f\. Moreover, it becomes clear that the normal TFR also 
holds for a more general form of the PDFs h(u,t ), that is, for a more general class of 
the positively valued stochastic processes r(u). 


Type A and B FFPEs: For these two types the fluctuation ratio in Laplace space is 
more complicated than for type C FFPE in Eq. 30. It is obtained with Eq. 24 as 


PaAW,s) = hAW/F^s) = f vo 
Pa,b(-W,s) / 4B i-H7/'„,») U,V 


(34) 


In contrast to Eq. 30, here the right hand side depends on the Laplace variable s. 
Consequently, one may expect an anomalous ratio R which is confirmed numerically in 
the overview of Fig. 2. The fluctuation ratios of type A (left column) and type B FFPEs 
(middle column) show a nonlinear increase as functions of W. For type B FFPEs there 
is a clear transition at the current maximum of the PDFs at W max = v a F Q t which is 
equal to t with u 0 = 1 and F 0 — 1 in Fig. 2. For W > W max the fluctuation ratio 
increases with time. In contrast, the TFR of type A FFPE increases faster in W than 
for type B. At the scale of this overview plot there is no transition point visible as for 
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type B FFPE. However, the qualitative time-dependence of the fluctuation ratio for 
type A FFPE is the opposite to type B FFPE: The ratio increases faster for smaller 
times. To gain further insight into this behavior, some asymptotic expansions of the 
TFR for type A and type B FFPEs are performed in the next section. 

3.3. Asymptotic expansions of the fluctuation ratio for type A and B FFPE 

In this subsection we analyze the asymptotic behavior of the work fluctuation ratio for 
type A and B FFPE. Differences between type A and type B simply correspond to the 
value of a which is 1 < a < 2 for the superdiffusive FFPE of type A and 0 < a < 1 
for the subdiffusive type B FFPE. Type C is not considered anymore, as the analytical 
calculation of Eq. 31 and the numerical analysis in Fig. 2 have delivered a normal 
fluctuation relation with a time-independent linear increase in the work W. 


Small W expansion: First, the behavior of the TFR for the work PDFs of the FFPEs 
is studied for small IE as a function of time. The logarithmic ratio of a continuously 
differentiable function p(z) can be expanded as Taylor series for positive z as 


log 


P(z) 

p(-z) 


2 dp(z) 
p(z = 0) dz 


z + 0(z 2 ). 

z =0 


(35) 


Inserting the approximate work PDF pfW, t) from Eq. 25 together with the 
transformation of Eq. 27 into Eq. 35 requires the calculation of the derivative of the Fox 
//-function. Using Eq. B.4 with r = 1, h = 1, c = —l/y/K a t a /F 0 , and d = v 0 t/y/K a t a 
allows us to calculate the linear term in the Taylor expansion of Eq. 35. With the 
assumption W/Fq < vpt and after some simplifications using the definition of the Fox 
//-function by the Mellin-Barnes integral in Eq. B.l one obtains the fluctuation ratio 
for small W as a quotient of two Fox //-functions: 


R(» 



ttI 0 
H 11 

v 0 1 

(1 - a/2, a/2) 


\J K a t a 

(i.i) 

W 

ttIO 

v 0 1 

(1 - a/2, a/2) 

Fo 

VK a t a 

(0.1) 



A (t) W. 


(36) 


The prefactor Aft) summarizes the time-dependence of the fluctuation ratio. Its 
numerical evaluation based on the Taylor series of Eq. B.3 is shown in Fig. 3A. In 
the superdiffusive case 1 < a < 2 (type A FFPE) the prefactor A(t) increases as a 
function of time, whereas in the subdiffusive case it decreases with time. The argument 
of the Fox //-functions z = v 0 t/\flKf t° in Eq. 36 scales ~ t 1 with 1 — a/2 > 0 for 
0 < a < 2. Thus the asymptotic expansion of these Fox //-functions can be used for 
t —> oo. In the long time limit the scaling function Aft) converges towards the following 
non-zero constant value: 


lim A (t) 

£—>•00 


2 | — 
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a/[a— 2 ) 
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1 /a 
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F 0 K\ f 0 


W for a = 1. 


(37) 
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Figure 2. Time dependence of the fluctuation ratio for type A FFPE (left column), 
type B FFPE (middle column) and type C FFPE (right column) for different values 
of a (rows) and times t = 1,2,4,8. Parameters were selected as K a = 1, v 0 = 1 
and A a vo = 1. Whereas a = 1 and all cases of a for type C FFPE show a normal 
fluctuation ratio with time-independent slope (in all of these cases the linear t = 8 curve 
hides the previous times t = 1, 2,4) all other sub-plots show a more complex time- and 
work-dependent fluctuation ratio: Anomalous non-Markovian dynamics and/or non- 
Gaussian behavior cause a complicated time-dependence and non-linear behavior of 
the work fluctuation ratio. 
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The corresponding values are shown as squares in Fig. 3A indicating the predicted 
asymptotic behavior. Fig. 3B shows the spatial behavior of the work fluctuation ratio 
for two subdiffusive examples a = 0.4 and 0.8 at different instants of time t (compare 
to small W values in the overview given in Fig. 2). The slope of the ratio decreases 
with increasing time and agrees well with the small W expansion given in Eq. 36. The 
superdiffusive case in Fig. 3C shows a reverse behavior as the small W ratio increases 
with time. As indicated in Fig. 3A it can also be negative as show in Fig. 3C for a = 1.6 
and t = 1,2. In the superdiffusive case, the small W expansion has a smaller region of 
agreement with the exact ratio. The more complex behavior is technically due to the 
two separating peaks of the PDF as illustrated in Fig. 1. 


Large W expansion: Finally, the behavior of the work fluctuation ratio is studied for 
large values of the work W. The overview given in Fig. 2 shows a different non-linear 
behavior for the subdiffusive and superdiffusive case. Assuming W/F 0 > Vot and large 
arguments of the Fox //-function for type A and type B FFPE in Eq. 25 allows us to 
use the asymptotic expansion of the corresponding Fox //-function in Eq. B.5. For large 
W one obtains the following relation: 


R(W,t)\ 


2 vpt 


W^-oo 


a \a/{2-a) 

2 


y/Kj c 


2 /( 2 -a) 


pp“/(2-a). 


(38) 


Thus the work fluctuation ratio scales as a power law with an exponent a/(2 — a). This 
exponent is between 0 and 1 for the subdiffusive type B FFPE. For superdiffusive type 
A FFPE it is larger than 1. This asymptotic power law behavior is shown in Fig. 4 for 
two examples. Continuous lines represent the result of Eq. 38 and agree for larger W 
values with the exact results denoted by circles. Eq. 38 additionally contains a time- 
dependent scaling factor that is proportional /( 2q_2 )/(“~ 2 ). This factor is positive for the 
subdiffusive type B FFPE and negative for type A FFPE. 


4. Summary and outlook 

In this work we studied three different types of FFPEs generating anomalous diffusion: 
a superdiffusive one (type A), a subdiffusive one (type B), and another one that exhibits 
a transition from sub- to superdiffusion under parameter variation (type C). Type A and 
type B break FDR1 while type C preserves it. Type A can be derived, under certain 
assumptions, from an overdamped Langevin equation with power law correlations of 
the velocity fluctuations, types B and C have been derived before in the literature from 
CTRW theory. Type C can also be obtained via subordination. We then calculated 
position PDFs for all models analytically and studied the shapes of all PDFs numerically 
under variation of the anomaly index as they evolve in time. Finally we checked the 
work TFR for all three models. Especially, we studied the time dependence of the ratio 
of the work fluctuations both for small and for large work by analytical asymptotic 
expansions in comparison to numerical evaluations. 
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work W work W 

Figure 3. A. Time dependent decay of the initial fluctuation ratio A (t) defined by 
Eq. 36 for small work W and different values of a corresponding to type A FFPE 
(1 < a < 2) and type B FFPE (0 < a < 1) with parameters K a = 1, Fo = 1, u 0 = 1 
and A a vo = 1. Circles show the direct calculation for small W from the ratio of PDFs 
as defined in Eq. 28 whereas lines result from the computation of the first term of the 
small W expansion of Eq. 35 and Eq. 36. Both calculations agree and A (t) converges 
towards the long time limit given by Eq. 37 as indicated by the squares. Whereas A (t) 
is time-independent for a = 1, it decrease or increases as a function of time for the 
subdiffusive (type B FFPE) and superdiffusive case (type A FFPE), respectively. B. 
The fluctuation ratio of work is shown for the subdiffusive case as a function of work 
and different time points as indicated. The slope decreases for increasing time. Thin 
black lines indicate the small work limit of Eq. 36. The obvious kink at W = 1 for 
t = 1 is due to the peak of the corresponding PDF in Fig. 1. C. The superdiffusive 
case shows a more complicated behavior: The small work slope increases with time. In 
addition, it also changes from negative to positive for small time in the a — 1.6 case. 
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Figure 4. Large W asymptotic of the work fluctuation ratio of type A and type B 
FFPEs. Continuous lines show the asymptotic large W result given by Eq. 38. Circles 
indicate the exact result from the direct computation of the work fluctuation ratio. A. 
Subdiffusive case for a = 0.8 corresponding to type B FFPE. B. Superdiffusive case 
for a = 1.2 as example for type A FFPE. 


We find that our type C model with FDR1 exhibits a conventional work TFR for 
all times, meaning the fluctuation ratio is constant in time and linear in the work. For a 
correlated Gaussian stochastic process it was shown that FDR1 implies the existence of 
a conventional TFR [43]. Our work generalises this result to an example of non-Gaussian 
PDFs generated by FFPE dynamics. It is interesting that the conventional TFR is still 
obeyed, despite the highly lion-trivial dynamics exhibited by both the position PDFs 
and the corresponding moments. The existence of the conventional TFR for this case 
is connected to the fact that only the equation for Type C describes a subordinated 
process, namely the one subordinated to Brownian motion with drift under random 
time transformation. An important open question is to which extent Fig. 1 in [43] 
summarising the interplay between FDR1, FDR2 and TFRs for correlated Gaussian 
stochastic processes in terms of necessary and sufficient conditions can be generalised 
to non-Gaussian processes. For our other two models type A and type B the position 
PDFs show also very subtle and non-trivial non-Gaussian shapes. However, in contrast 
to type C they are characterised by a highly non-trivial fluctuation ratio: For type 
A the latter decreases with time, for type B it increases. Similar results have been 
obtained for the work TFR of strongly correlated Gaussian stochastic processes without 
FDR1 [42, 43]. On top of this, for both types of FFPEs the fluctuation ratio yields 
different long time limits depending on whether the work is small or large: For small 
work the fluctuation ratio converges to linearity in the work with constant prefactors, 
which reminds of the conventional TFR; however, here the slopes depend on the anomaly 
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index of the dynamics. For large work the fluctuation ratio remains nonlinear in the 
work, with convex and concave shapes for type A and type B, respectively. 

Our work was motivated by experiments on cell migration [50], where data were 
successfully fitted by solutions of a fractional Klein-Kramers equation [48]. Several 
generalisations of such a Klein-Kramers equation have been proposed to describe 
processes under external fields [48, 49, 51], which in turn yield FFPEs for the position 
only, similar to the ones studied in our paper, as special cases [21, 52, 12], We thus 
believe that our present work might have important applications to understand cell 
migration in nonequilibrium situations such as under chemical gradients; see [44] for 
first results. More generally, our theory might have applications to understand glassy 
nonequilibrium dynamics: In computer simulations of a number of glassy systems 
violations of conventional TFRs have been observed featuring fluctuation ratios that 
are nonlinear in the work with time-dependent prefactors [40, 41], 

Apart from such experimental applications, our first approach for deriving a FFPE 
pioneered by Balescu [57, 58] deserves to be studied in more detail. For example, it 
would be interesting to derive a superdiffusive FFPE from it that preserves FDR1, and 
to check again the TFR. On a broader scale it would be important to generalise our 
approach by considering more general observables, ideally dissipation functions [1] or 
related functionals defined within stochastic thermodynamics [7]. More general force 
fields than simply constant forces [42] and other types of FRs could be tested as well. 
Such theoretical studies may pave the way to identify different classes of anomalous FRs 
characterized by specific functional forms, generalized FDRs associated with them, and 
to explore the physical significance of these results. Last not least the quality of the 
Galilean invariant approximate solution Eq. 25 [20, 12] of the FFPEs 8,10 needs to be 
investigated in detail. 
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Appendix A. Pseudo-Liouville approach 

Following the so-called pseudo-Liouville hybrid approach of Balescu [57, 58] allows us to 
relate the dynamics of a particle defined by a Langevin equation to the corresponding 
PDF of the stochastic process. We start from the Langevin equation for the position 
x(t ) of a particle 

^ = «. + »(*). (A.1) 

where v(t) is a correlated stochastic process with zero mean (v(t)) = 0 and a given 
correlation function (v(t)v(t')) = T(t — P), where the average is performed over the 
stochastic process v(t). Vo denotes a constant external force. The stochastic function 
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F(x,t) 

F(x, t) = 5{x — x(t)) (A.2) 

represents the exact density of the process. Derivation of Eq. A.2 with respect to time 
and the usage of the Langevin equation Eq. A.l delivers the continuity equation for the 
exact density F(x,t): 


dF(x,t) d ... dx(t) 

-- x(t)) —f 1 

dt dx dt 


dF(x,t ) , dF(x,t) .. . 

1 ; +[uo + u(t)] l ’ =0. (A.3) 


dt dx 

Now, the exact density F(x,t) is decomposed into an averaged part f(x,t ) and 
fluctuations 5f(x,t ) 


F(x,t) = with f(x,t) = (F(x,t)). 


(A.4) 


ft is the further aim of this appendix to calculate the PDF f(x,t ) for the stochastic 
process defined by the Langevin equation Eq. A.l for given correlations of v(t). 
Averaging of the exact density in Eq. A.3 leads to 


{m +v i) nx ' t) = -rf v(t) 6 f(x F- 


(A.5) 


(A.6) 


Subtraction of Eq. A.5 from Eq. A.3 results in 

^ + [ y o + f( x ^) _ 

Eq. A.6 can be solved with the method of characteristics 

t 

= J dt'(v(t') f(x - A(t,t'),t') - (v(t') Sf(x - (A.7) 


with the definition A (t,t') = Vo(t — t') — J dt\ v(ti). Inserting Eq. A.7 into Eq. A.5 

t' 

delivers the final equation for the PDF f(x,t): 


d d \ d 2 k 

^ + V °dx) ^ X,t ^ = dx 2 I dt '( v ^ v ^ f( x ~ A (M')>0)- 


dt 


(A.8) 


This is an exact relation for f(x,t) that is generally non-local in space and non-local in 
time, i.e. non-Markovian. Applications and approximations of this relation are studied 
in section 2.1. 


Appendix B. Definiton and properties of Fox 77-functions 

The Fox 77-function is defined as inverse Mellin transform of the function y(s) [12, 71] 


h^{z) = h;] 
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(bj, Bj)j = i q 


= 2^ii x(s)z ‘ is 
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11 r 0j - B i s ) 11 r (l - % + Ajs) 

X(s) = -j=±-p-, (B.2) 

n r(i - bj + Bjs ) n r(aj - Aj-s) 

j=m+l j‘=n+l 


0<n<p, l<m<q } (aj, bj) e C., and (Aj, Bj) e M + . Empty products in Eq. B.2 are 
taken as one. 

A series expansion allows the numerical calculation of Fox //-functions. The 
following form for a special Fox //-function is used: 
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(B.3) 


Summation in this work is performed numerically with multiple-precision arithmetic. 

The derivation of the Fox //-function is required to calculate the fluctuation ratio 
for the FFPEs of type A and B. This can be performed using the following relation [72]: 
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(B.4) 


For large arguments the Fox //-functions of type H^(z) decay as stretched exponential 
functions. The asymptotics of the PDF in Eq. 25 is given for large z by [73, 72] 


f A,B (■*) t) ~ 


1 — a 

Z 2 -« 

a/47t(2 — a)K a t a 



exp 


for z —y oo and z — \x — v 0 t\/\/K a t a . 



(B.5) 
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